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Abstract 

Modal Identification results are presented 
that were obtaired from recent flight flutter tests 
of a drone vehicle with a research wing (OAST ARW-1 
for Drones for Aerodynamic and Structural Testing, 
Aeroelastic Research Wing-1). This vehicle is 
equipped with an active flutter suppression system 
(FSS). Frequency and damping of~severa1 modes 
are determined by a time domain modal analysis of 
the impulse response function obtained by Fourier 
transformations of data from fast swept sine wave 
excitation by the FSS control surfaces on the wing. 
Flutter points are determined for two different 
altitudes with the FSS off. Data are given for 
near the flutter boundary with the FSS on. 


I. Introduction 

One application of active controls technology 
to aircraft design is to suppress flutter. The 
potential benefits have been well demonstrated in 
wind tunnel tests and by analysis (e.g.. Refs. 1 
and 2). A flight test program has been initiated 
by NASA to demonstrate active control technology 
including flutter suppression systems (FSS) in 
flight using a remotely piloted vehicle or drone 
aircraft. Free flioht data can thus be acquired 
on an aircraft wi' educed risk. This program is 
called OAST (Dron* .or Aerodynamic and itructural 
Testing) and is a cooperative effort by the NASA 
Langley and Dryden Research Centers with contractual 
support front the Boeing Military Airplane Company, 
Wichita, Kansas. An overview of the OAST program 
is given in Ref. 3. 

Three flights have been flown with the first 
aeroelastic research wing (ARW-1) from October 1979 
to June 1980. Due to several system problems 
including a hydraulic pump failure and difficulties 
with the telemetry systems, the first flight 
(October 1979) was essentially a systems development 
test. During the second flight (March 1980) suf- 
ficient structural dynamic response measurements 
were made with the flutter suppression system off 
(FSS OFF) to obtain a good estimate of the flutter 
Mach number for cne altitude. The third flight 
(June 1980) included flutter testing with the FSS 
on and off. During this flight, flutter was en- 
countered with the FSS on p-imarily as a result of 
an erroneous setting of the FSS gain to one-half the 
desired value. An ove-all description of these 
testssii Oil -1 in’ testi g procedures, a description 
of the flutter incident, and some of the results 
of these flights are presented in Ref. 4. 

The flight time of the OAST vehicle is limited 
to about 20 minutes. T.his short flight time and 
‘he additional complication of testing with the FSS 


both on and off required the development of a rapid 
testing technique. Fast-swept sine wave excitation 
by the FSS controls is used at each flight condition 
to yield data to define the dynamic characteristics 
of the vehicle. These dynamic data were analyzed 
on-line in near real time to detennine the damping 
and frequency of the principal mode (Ref. 4). Sub- 
sequent data analysis has been performed after each 
flight to define additional modes. This paper 
presents the results of these post-test or off-line 
analyses of the fast swept sine data using fast 
Fourier transform techniques along with a time- 
domain modal analysis^ to extract the modal fre- 
quencies ana dampings. A brief description of the 
overall test procedure is also piesented and the 
modal analysis procedure is described. Flutter cal- 
culations for the DAST vehicle are compared in 
Ref. 6 with the experimental results of this paper. 


II. Nomencla t ure 

a^, coefficient in curve fit. the offset or static 
value (Eq. (1 )) 

a. coefficient or kth cosine term in curve fit 

(Eq. (D) 

b. coefficient of kth sine term in curve fit 

(Eq. (D) 

a^ acceleration, g's 

E mean-squared error (Eq. (2)) 

E expected value (Eq. (3)) 

F Fourier transform 

f frequency, Hz 

h altitude, km (ft) 

i data point index, 1 to N 

j parameter index 

K number of modes in curve fit 

k modal index, 1 to K 

M Mach number 

N number of data points in digitized time 
history 

output error covariance matrix 

3 parameter sensitivity matrix 

t time, seconds 

V velocity 
flutter velocity 

Y curve-fitting expression (Eq. (1)) 

Ap increment in parameter vector 

; fraction of critical damping 

S control deflection, degrees 

n damping coefficient (Eq. (1)) 

i. frequency, rad/sec 

Su bscripts 

1 left wing 

r right wing 

c command signal 
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111. Description of Vehicle and Test s 

The DASl ARW-1 dircraft. shown In Fig. 1, 

IS d modi tied ii rebee II drone. It has an overall 
length of 8.6 meters (28.3 ft). The ARW-1 research 
wing was designed for cruise near M • 0.98 and 
has a span of 4.34 meters (14.3 ft), a supercritical 
airfoil section, and an aspect ratio of 6.8. The 
wing was designed to flutter within the flight 
envelope in order to assess the effectiveness of 
the onboard active control system used to suppress 
flutter. 

The OAST operational concept is shown schemati- 
cally in Fig. 2. The drone is carried aloft under- 
neath the wing of a B-52, air launched, and flown 
remotely by a pilot in a ground cockpit. After a 
flight of 20 to 30 minutes, the drone is recovet'ed 
by mid air capture and retrieval by helicopter 
(Refs. 3 and 4 gives further description of the 
operational aspects). 

The OAST vehicle is extensively instrumented 
for measurement of rigid body type parameters, 
loads by calibrated wing strain gage bridges, 
wing steady state surface pressures, and the flutter 
system parameters. The data are digitized onboard 
and transmitted by two telemetry systems in pulse 
code modulated form to the ground for on-line use 
and recording. The sampling rate of the various 
channels varies from 20 to 500 samples per second. 
Anti-aliasing filters were incorporated onboard 
for each channel, ihe high sample rate signals 
(500 sps) are used for the flutter system parameters 
and include the FSS accelerometers, control position 
potentiometers, and the control system excitation 
signal. These data are recorded and later converted 
to a digital tape for computational analyses. 

A sketch of the flutter suppression system for 
the DAST ARW-1 is shown in Fig. 3. The FSS includes 
wing tip and inboard accelerometers, a hydraulic 
control system with fast response actuators, wing 
aileron control surfaces, and onboard electronics. 

The FSS ailerons are used for flutter suooression 
and dynamic excitation but are not used for flight 
contix)! of the aircraft. The system is designed 
to suppress both synmetne and antisynmetric flutter 
and can be remotely engaged or disengaged by ground 
command in the structural analysis facility (SAr) 
where the flutter experimenters are located (Fig. 2). 
Since the flutter frequencies are in the 10 to 
20 Hz range, the FSo control system is designed to 
respond at high frequencies. The lowest natural 
frequency for the control surface, actuator, and 
hydraulic system combination is near 50 Hz. This 
fast response aileron control system is also used 
for structural dynamic excitation. Control surface 
pulses (a single cycle of a 20 Hz sine wave) are 
used for qualitative cn-line monitoring of aero- 
elastic stability particularly during accelerating 
flight. For a more complete evaluation of stability, 
a fast swept sine wave input by the controls is 
used. The swept sine signal runs from 10 to 40 Hz 
in 7 seconds and is tapered near the ends of tne 
sweep to minimize transients as suggested in Ref. 7. 
The sweeps of 7 seconds duration are considerably 
shorter than normally used in the flutter testing 
of aircraft. The short flight duration in combina- 
tion with the need for synmetric, antisymmetric. 

FSS on, and FSS off testing led to the selection 
of the short duration of the sweeps. The relative 
high flutter frequency of 10 to 20 Hz also suggests 


that the sweep duration could be shorter than 
usually used for lower frequency modes. The sweeps 
and poises are generated onboard the aircraft and 
are commanded from the ground by an operator in the 
SAF. Both symmetric and antisynmetric inputs can 
be generated. For test points of nominally constant 
flight conditions, a series of synmetric and anti- 
symmetric sweeps and pulses are used. The data from 
the sweeps are analyzed on-line in near real time 
to obtain an indication of flutter mode stability 
(Ref. 4). For points below the predicted flutter 
boundary, the excitation is performed with the FSS 
both on and off. After the FSS off flutter boundary 
is approached or exceeded, the FSS is left on and 
the testing is performed only for the system on 
condi tion. 


Pata A na lysi s 

Overal 1 Procedure 

The block diagram of the post-test analysis of 
thf ta with fast swept sine wave inputs is 
presented in Fig. 4. The objective is to obtain 
gooo estimates of the frequency and damping of the 
modes of motion from the measured input and output 
signals. The procedure used here involves taking 
Fourier transforms of these signals numerically, 
forming the frequency response functions by dividing 
the transforms, and after a filtering operation, 
taking the inverse Fourier transform to obtain the 
impulse response function. The impulse response 
function contains the free vibration characteristics 
of the modes that respond to the excitation. A 
time domain modal analysis procedure suitable for 
analysis of free vibration motion (Ref. 5) is 
used to determine the frequency and damping of the 
principal modes. There are several alternative 
methods available for modal analysis (e.g. , Refs, 5, 
10, and 11) depending upon the type of data (fre- 
quency domain, etc.) available and the desired 
form of the results The method used is one that 
was designed for this type of data (Ref. 5) and 
is conceptually simple. The data analysis procedure 
(Fig. 4) was implemented on the CDC CYBER 175 com- 
puters at Langley Research Center. The steps in 
this procedure are discussed in further detail in 
the following paragraphs. 

The first step in the procedure is to sum the 
output of the left and right FSS accelerometers 
to minimize any extraneous antisymmetric motions 
occurring during symmetric excitation (Fig. 4) and 
to difference the accelerometer signals for anti- 
syimietric excitation (as recommended in Ref. 8). 

The input signal for analyzing the data for FSS off 
cases consists of the control deflection signals 
also surnned or differenced depending on whether the 
motion is synmetric or antisynmetric. With the FSS 
on, the control deflection includes the sum of 
deflections in response to the excitation conmand 
and the deflections for flutter suppression. Use 
of the control deflection as the reference input 
signal gives the system off dynamics. For analysis 
with the FSS on, the sweep excitation conmand signal 
is used as the reference input signal (Fig. 4). 

The discrete finite Fourier transforms of the 
input signal and the output signal (accelerometers) 
are calculated and the frequency response function 
is formed as the ratio of the two transforms As 
previously indicated, the swept sine wave input ran 



from 10 to 40 Hz. Although the frequency response 
function Is calculated from 0 to 250 Hz (the Nyquist 
frequency), there is essentially no input power and 
the frequency response function is meaningless 
outside the 10 to 40 Hz frequency range. Hence, 
the frequency response function is filtered by 
multiplying by a window function in the frequency 
domain. The window function used is unity from 
10 to 37.5 Hz, and is tapered with a 1-cos^ function 
from 7.0 to 10 Hz and from 37.5 to 42.5 Hz, and 
is zero at all other frequencies. The impulse 
response function is generated by taking the inverse 
Fourier transform of the filtered frequency response 
function. The filter function eliminates both 
spurious frequency response function estimates where 
the input power is nearly zero and any real noise or 
data in the frequency ranges where the filter func- 
tion is zero. One consequence of band-limiting the 
frequency response function in this manner is a 
smearing of the impulse response function. The 
resulting impulse response function is a circular 
convolution of the transform of the filter function 
and the true impulse response function. This effect 
primarily produces a smearing phase error near the 
t = 0 portion of the impulse response. The l-cos^ 
tapering decreases the side lobes of the time 
domain window thus reducing the smearing effect. 

This tapering was chosen to be a compromise between 
noise elimination and windowing effects. The 
windowing effect is also partially alleviated by 
analyzing the impulse response starting at t = 0.05 
(instead of t - 0) and by allowing additional terms 
in the complex exponential time domain fitting 
technique. Generally a moderate amplitude term 
near the lower cutoff frequency (8 to 9 Hz) appears 
in the results and a small amplitude term near 40 Hz 
is occasionally obtained. 

A similar data analysis procedure was used on- 
line for analyzing the fast swept sine wave data 
with a minicomputer (Ref. 4). The on-line modal 
analysis assumed a single degree of freedom model 
due to the time constraints of the f^light, whereas 
the procedure used here takes into account all the 
principal modes in the response. 


Inspection of Eq. (1) shows that if nk and 
uk and the value for K are known, it is possible 
to compute ag, ak, and bk by solving a linear 
least-squares prob em. The parameters nk <^k 

enter the minimization problem in a nonlinear manner 
and must be determined by some type of search algo- 
rithm. Although this is a standard nonlinear 
optimization problem for which several methods are 
available, for simplicity a direct search technique 
is used to search the coordinate space (nk> <ak) 
until a minimum is neared. A linearized tecimique 
is then used to accelerate convergence to the mini- 
mum. The method requires initial estimates for 
nk and uk. for a single mode the inital 
estimates can be nearly arbitrary. However, for a 
multiple-mode case, the computer time can be signif- 
icantly reduced by choosing good initial estimates. 
The following procedure has been found to be a 
reasonable way of getting the initial estimates for 
a multiple-mode case: 

(a) Generate a one-mode solution using arbi- 
trary initial estimates. 

(b) Compute the difference between the one- 
mode solution and the input data, that is, the 
output error. Then generate a ono-mode fit to the 
error. 

(c) Use the ok and ^k values from steps 
(a) and (b) as the initial estimates for the two- 
mode solution. 

(d) For higher modes, steps (b) and (c) arc 
repeated using the difference between the current 
multi-mode solution and the original data to 
estimate the next higher mode. This procedure is 
continued until enough terms (K) are used to make 
the output errors small. This method has been 
Implemented in a research code on the CDC CYBER 175. 
A typical run of a four mode fit to 256 points in 
the time history takes about 50 to 60 seconds of 
CPU time including plotting. 

Estimates of Standard Deviations of the Parameters 


Modal Analysis Technique 


The time domain modal analysis procedure ana- 
lyzes a single channel of data. Given a free decay 
record such as the impulse response function which 
contains the response of one or more vibration modes 
in the form of a digitized time history, the problem 
is to determine the damping, frequency, amplitude, 
and phase of each mode. A least-squares curve fit 
of the data is made with complex exponential func- 
tions (or damped sine and cosine waves) in the form 


Y(t) = Oq + 



cos U|^t + b|^ sin uj^t) 


( 1 ) 


The parameters an, av, b^, nki and uit are 
determined such that the squared-error difference 
between the expression Y(ti) and measured time 
history a^ is minimized. The error is given by 

N 

E=J][Y{t.) -a^^]2 (2) 

1 = 1 


One difficulty with parameter estimatin'' from 
dynamic data is determining the quality of the 
estimate. Usually repeatability, scatter, etc., 
give the experimenter some indication of the quality 
of results. Another approach is to estimate the 
standard deviations of the determined parameters 
(Ref. 5). The standard deviations of the estimated 
parameters, or uncertainty levels, can be determined 
from maximum-likelihood theory (Ref. 11, e.g.) 
based on the Cramer-Rao bound. This type of esti- 
mate has provided some useful results in the field 
of stability and control (e.g.. Ref. 9). Assuming 
only measurement noise that is Gaussian and white, 
the expected variance of the parameter vector is 


E[a5 ap'] = 


N r 


i = lL 


S^(t.)R,'^ §(t^) 


(3) 


where S is the parmeter sensitivity matrix, R] 
is the output error covariance matrix, (here a 
constant) and T denotes matrix transpose. The 
parameter vector p is made up of ag. a^, b^, 
nk, and u>k. The parameter sensitivities are 
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3Y(t.) 

given by Sj(t^) = -jp — and can be calculated 

by directly differentiating Eg. (1). The variance 

is a nonnal output parcmeter of the program. 
Thus, for a single channel of data, as considered 
here, these parametc" us ertainty levels can be 
readily calculated arter the curve-fitting process 
is completed. 

One difficulty with the estimates of the stan- 
dard deviations in the pasc has been that they gave 
unrealistically low values. This is discussed in 
Ref. 9 and is attributed to the fact that Eg. (3) 
assumes a uniform error spectrum, whereas in prac- 
tice the error spectrum is band limited. In the 
present case the errors are largely a result of 
background modal response. Although some allowance 
for the actual spectrian of the error can be made, 
it is suggested (Ref. 9) that multiplying the result 
calculated from (3) by 5 to 10 gives realistic 
results. A factor of 10 is used in this report. 
Experience has shown t‘'at these estimates are then 
a good gualitative guide to the validity of the 
estimated parameters. 

Saimjle Case 

An example of a time history recording from the 
OAST flight test is shown in Fio. 5. The resoonse 
is from FSS jff, symmetric excitation at M = 0.74 
and h = 4.6 m (15.000 ft). The resulting impulse 
response function and a four mode fit is shown in 
Fig. 6, and the modal content of the fit is shown 
in Fig. 7. The low freguency mode labeled mode 4 
is a result of the window function as previously 
discussed and is not a structural vibration mode. 

The three modes are identifiable by comparing with 
the ground vibration modal freguencies of Table 1. 
The first synmetric wind bending is mode 1, a 
synmetric wing- fuselage mode is mode 2, and mode 3 
is the first symmetric torsion mode. The noise level 
both in the measured data (Fig. 5) and in the im- 
pulse response (Fig. 6) is very low. 


V. Results and Discussion 

Data from two flights showing the approach to 
the flutter boundary at two diffei-ent altitudes 
will be discussed. The altitude and Nach number for 
the test points are shown in Fig. 8. The approach 
to near flutter for flight 2 was at a nominal 
altitude of 7.6 km (25,000 ft) with the FSS off. 

Both syimietric and antisymmetric data are presented 
for flight 2. For flight 3 the approach to flutter 
was for an altitude of approximately 4.6 km 
(15,000 ft). The data were taken with the FSS on 
and off at the lower Mach numbers, but only wich 
the FSS on at the higher Mach numbers (Fig. 8). 

For flight 3, only data for the syitiretric modes 
are presented since they were critical for flutter 
with the FSS both on and off. 

Tests at 7.6 km (25,000 ft) Altitude (Flight 2) 

The estimated damping and freguency for the 
first bending and torsion modes for both synmetric 
and antisynmetric motion are shown in Fig. 9. The 
damping for the bending mode approaches zero and 
flutter is indicated near M : 0.92 for both 
symmetric and antisynmetric bending modes. There 


is little scatter apparent in the damping or fre- 
guency of the flutter mode and the flutter condition 
is approached with a moderate rate of decr-Ms.? in 
damping with Mach number. The results for the 
torsion mode shows some scatter particularly in the 
damping values. The dampin' of the torsion modr 
increases and the freguency decreases as flutter is 
approached. These trends are typical for a hending- 
torsion type flutter as the flutter boundary is 
approached. 

The data from Fig. 9 are repeated in Fig. 10 
and the calculated standard deviation estimates are 
shown as the vertical bars through the symbols. 

The results indicate that the bending (flutter) mode 
is well estimated but that the torsional mode is not 
as well defined. This result is apparently a 
result of the higher damping and smaller amplitut,.' 
of the torsion mode. The standard deviation esti- 
mates are consistent with the level of scatter and 
are smaller for the freguency than for damping. The 
result that freguencies are estimated more accurately 
than damping is generally the case for ‘■‘lis type of 
data. 

Four modes were used in fitting the impulse 
response function for these results. There is a 
fuselage first bending mode in the frequency range 
(Table 1), but it could not be reliably extracted 
from the FSS accelerometer data. Further definition 
might be possible from other data channels such as 
the inboard accelerometers. 

The results for the flutter mode indicate a 
low level of scatter and also correspond guite well 
to the on-line results (Ref. 5). This low level 
of scatter is likely a result of the low noise level 
in the data. 

Test s at 4. 6 km (15,000 ft) Altitude (Flight 3) 

FSS Off . The symmetric damping and freguency 
estimates determined for the FSS off are shown ii 
Fig. 11. Three modes are indicated which correspond 
to the first wing bending, first fuselage bending, 
and first torsion. Flutter near M = 0.80 is 
indicated in the mode associated with first syin- 
metric wing bending. For this altitude there is j 
large interaction of the bending and fuselage modes 
and their frequencies apparently cross as Mach 
number is increased (Fig. 11). Dynamic dat' with 
response in two or more modes of nearly the oame 
frequency can be troublesome for obtaining good 
estimates of the modal parameters. The trends shown 
by the lines in Fig. 11 are based primarily on 
continuity of the damping values and are possibly 
subject to some interpretation. These results were 
obtained using a four-mode approximation and the 
FSS accelerometer data. 

The standard deviation estimates are shown in 
Fig. 1? for each of the three modes. The parameters 
are generally well estimated except for Mach numbers 
near 0.77 where the freguencies are nearly coinci- 
dent. The fuselage modal participation was apparent- 
ly sufficiently large at this altitude that the 
damping and frequency could be estimated. The 
damping for the fuselage mode is low and the stan- 
dard deviation would indicate the possibility of 
instability. This mode was certainly stable for the 
FSS off condition at M = 0.70 and M = 0.74, and 
the actual value of damping is more likely that of 
the symbols. 



For the two Mach numbers for which measurements 
were made with the FSS both off and on, M = 0.7 
and 0.74, the results for the FSS off are in very 
good agreement (Fig. 11). This good agreement also 
probably results from the low noise level in the 
measured data. 

FSS On . The syimetric damping and frequency 
values determined for two of the principal modes 
with the FSS on are shown in Fig. 13. These two 
modes correspond to wing first bending and first 
torsion. The frequency of the bending mode is 
significantly raised by the FSS from 12 to 15 Hz to 
near 20 Hz. The damping is also significantly 
increased over that measured with the FSS off at 
the lower Mach numbers, but flutter would be pro- 
jected to occur near M = 0.80 for the gain 
setting of the feedback system actually used for 
this flight. The aircraft was accelerating when the 
onset of sustained oscillations was obtained near 
M = 0.817. It had further accelerated to M = 0.827 
when the divergent oscillations led to break up of 
the right wing. The on-line da.mping estimates did 
not indicate the rapid onset of flutter as shown in 
Fig. 13. Further discussion of the flutter incident 
in given in Ref. 4 and comparisons with analytical 
results are given in Ref. 6. 

The estimates of the standard deviations of the 
damping and frequency values are shown in Fig. 14. 
The values for the bending mode are well determined 
except for the point at M = 0.70 which has a large 
standard deviation, possibly as a result of the 
higher damping value. The results for these two 
modes were determined from an eight mode fit to the 
data, and are the two modes that were consistently 
determined from the response. 


VI. Concluding Remarks 

The results of off-line analysis of data from 
two flights of the DAST ARW-1 have been presented. 
The dynamic rharacte'^istics of several iuodes have 
been obtained using the data analysis technique. 

This technique consisted of fast swept sine excita- 
tion by the aileron controls ^nd time-domain anal- 
ysis of the impulse response function derived by 
numerical Fourier transform techniques. The 
approach to the flutter boundary is clearly defined 
by the results of this technique. A significant 
contributing factor in the success of this technique 
is the low level of noise in the measured data. 

The calculated standard deviation of the parameters 
based on the Cramer-Rao bound gives a useful qual- 
itative indication of the quality of the frequency 
and damping estimates. 
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Tjble I Nomul mode frequencies predicted by 
NASTRAN analysis and measured during ground 
vibration test 
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Fig. 1 HAST ARW-1 vehicle In flight 


Fig. 4 Block diagram of post-test data analysis 
procedure. 
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Fig. 5 Sample input control deflection and accel- 
erometer responses for FSS off, M » 0.74, and 
h = 4.6 km (15,000 ft). 


Fig. 7 Components of sample impulse response 
function for FSS off, H * 0.74, h • 4.6 km 
(15,000 ft). 
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Fig. 6 Sample impulse response function and fit Fig. S Altitude and H.Kh number for test points, 

for FSS off, M » 0.74, h = 4.6 km (15,000 ft). 
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Damping and frequency estimates for FSS 
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(b) Anti syimetric 
Fig. 9 Concluded. 


Fig. 10 Standard deviations of frequency and damp- 
ing estimates for FSS off, h = 7.6 km (25,000 ft). 
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Fig. 10 Concluded. 
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Fig. 11 Synmetric modal frequency and damping 
estimates for FSS off, h “ 4.6 km (15,000 ft). 
Flagged symbols are FSS off results calculated from 
FSS on data. 
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Fig. 12 Standard deviations of frequency and 
damping estimates for FSS off, h = 4.6 km 
(15,000 ft). 
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Fig. 14 Standard deviations of frequency and 
damping estimates for FSS on. h = 4.6 km 
(15.000 ft). 
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estimates for FSS on, h * 4.fi km (b.OOO ft). 
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